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Abstract 



Bremsstrahlung radiation in the deuteron - proton collision 



Despite the fact that Bremsstrahlung radiation has been observed many years ago, 
it is still the subject of interest of many theoretical and experimental groups. Due 
to the high sensitivity of the AW — > NNy reaction to the nucleon - nucleon potential, 
Bremsstrahlung radiation is used as a tool to investigate details of the nucleon - nucleon 
interaction. Such investigations can be performed at the cooler synchrotron COSY in 
the Research Centre Jiilich, by dint of the COSY-1 1 detection system. 
For the first time at the COSY-1 1 experiment signals from y- quanta were observed in 
the time - of - flight distribution of neutral particles measured with the neutral particle 
detector. 

In this thesis the results of the identification of Bremsstrahlung radiation emitted via 
the dp — > dpy reaction in data taken with a proton target and a deuteron beam are 
presented and discussed. 

The time resolution of the neutral particle detector and its timing calibration are crucial 
for the identification of the dp — > dpy reaction. Therefore, methods of determining the 
relative timing between individual modules - constituting the neutron detector - and of 
the general time offset with respect to the other detector components are described. 
Furthermore the accuracy of the momentum determination of the registered neutron 
which defines the precision of the event reconstruction was extracted from the data. 
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I. Introduction 



In collisions between nucleons electromagnetic radiation can be emitted due to the 
rapid change of the nucleon velocity. This radiation is referred to as bremsstrahlung 
radiation. Although it has been observed many years ago, it is still the subject of interest 
of many experimental and theoretical groups [1, 2, 3, 4]. Due to the high sensitivity of 
the AW — > AWy reaction to the nucleon-nucleon potential, bremsstrahlung radiation is 
used as a tool to investigate details of the NN interaction and to discriminate between 
various potential models [1,2]. 

Such investigations can be performed at the cooler synchrotron COSY in the Re- 
search Center Jiilich, by means of the facilities installed there like the internal COSY-1 1 
detection setup. Although the COSY-1 1 facility is designed for threshold meson pro- 
duction covering a small solid angle in the laboratory system the special configuration 
allows also the study of reaction channels at high excess energies. The main aim of 
this thesis is an identification of bremsstrahlung radiation in a data sample taken by 
the COSY-1 1 collaboration during a measurement in January 2003 [5, 6]. The mea- 
surement was carried out at a hydrogen cluster target [7] using a deuteron beam with 
a momentum of 3.204 GeV/c. The registration of gamma quanta became possible be- 
cause the COSY-1 1 detection system has been extended by a neutral-particle-detector. 
This detector enables not only to measure the bremsstrahlung radiation created in the 
collision of nucleons but also opens wide possibilities to investigate the isospin de- 
pendence of the meson production in the hadronic interactions [8]. For example at 
present the COSY-1 1 facility permits to investigate T]-nieson production in proton- 
proton and proton-neutron collisions [9, 10]. Neutrons and gamma quanta registered 
in the neutral-particle-detector are identified via the time-of-flight on the 7 m distance 
between the target and the detector. In case of a neutron the time-of-flight and the hit 
position enables to determine its four-momentum vector. Since the time resolution of 
the detector and its timing calibration are crucial for the identification of the studied 
reactions, the data analysis of the dp — > dpy reaction will be preceded by the detailed 
presentation of the method used for the time calibration of this detector. 

The present thesis is divided into six chapters. The second chapter — following 
the introduction — describes briefly the motivations to investigate the AW — > AWy 
process, in particular, in view of the study of the proton-T] interaction by the COSY-1 1 
group. 

Description of the detection system, with a special emphasis on the structure of the 
neutron detector and the basis of its functioning will be presented in the third chapter. 

The fourth chapter is devoted to the time calibration of the neutral-particle-counter. 
The method of determining the relative timing between the individual modules and the 
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general time offset with respect to the other detectors will be depicted. This chapter de- 
scribes also the determination of the momentum resolution of the registered neutrons. 

The results of the identification of the dp — > dpy reaction in data taken in January 
2003 are presented in the fifth chapter, where also a detailed description of the data 
analysis is included. 

The sixth chapter comprises summary and perspectives, in particular a possibility 
to study the pentaquark state © + at the COSY-11 facility is discussed. 



2. Motivation 



For many years experimental and theoretical studies have been devoted to distinguish 
among various potential models and estimate of the off-shell amplitudes. The idea 
of using nucleon-nucleon bremsstrahlung as a tool for investigating this problem has 
attracted attention since calculations of the cross sections of the AW — > NNy reaction 
are highly sensitive to the AW potential. 

We would like to study this process also in view of another important issue related 
to the T]-proton interaction. 




Figure 2.1: Differential cross section for the p p — > p p r\ reaction as a function of the invariant 
mass of proton-proton system. Data are from ref. [11] and the line corresponds to the model of 
reference [12]. Picture adapted from [11]. 

Figure 2.1 shows a differential cross section for the pp — > ppr\ reaction as a func- 
tion of the invariant mass of the proton-proton system. The measurement has been 
performed at an excess energy of Q = 15.5 MeV. The theoretical description does not 
agree with data in the whole range of invariant mass of the proton-proton system. It 
has been pointed out that the difference between the data and the theoretical predic- 
tions originates from proton-r] interaction in the final state. In order to prove that this 
conclusion does not depend on the applied model it would be desirable to compare 
this distribution with the one obtained for the ppX system, where X does not interact 
strongly with protons. Therefore best suited for this purpose would be the ppy system. 
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Let us assume, that the differential cross section for the pp — > ppy reaction as a func- 
tion of proton-gamma invariant mass are determined. If the theoretical model, used for 
the calculations presented by the solid line in figure 2.1 were in perfect agreement with 
those data, it would corraborate the assertion that the disagreement between data and 
theory in the ppr\ case is due to the p-T] interaction. 

Studies of the pp-bremsstrahlung by means of the COSY-1 1 facility haven't been per- 
formed, but production of bremsstrahlung radiation in proton-proton collisions was 
investigated at the COSY-TOF facility [3]. Data have been taken at a proton beam mo- 
mentum of 797 MeV/c using a wide angle spectrometer. At the COSY-1 1 facility as 
a first step of investigations of the bremsstrahlung radiation the dp — > dpy reaction is 
studied, and the main point of the present work is an identification of dpy events in data 
taken during the experiment devoted to the measurement of the T]-nieson production in 
deuteron-proton collisions. 



n 1 
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Figure 2.2: Total cross section for the dp — > dpy reaction as a function of the excess energy. 
Original data are taken from [4]. 

A similar experiment — measurement of the dp — > dpy reaction — has been 
performed by the WASA/PROMICE collaboration at the storage ring CELSIUS, with 
deuteron beam energies between 437 MeV and 559 MeV. In these studies angular and 
spectral distributions are divided into two groups that can be attributed to a quasifree 
np — > npy process, and to a gamma production process with all three nucleons in- 
volved, viz: dp — > dpy. The results of these investigations — the total cross section 
for the dp — > dpy reaction — are given in figure 2.2. 

In this work the feasibility of the bremsstrahlung measurement at the COSY-1 1 
detection setup will be proven, and as soon as the analysis is finished the total cross 
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section data-base for the dp — > dpj reaction (fig. 2.2) will be extended by one point 
at Q = 557 MeV. 
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3. Experimental setup 



3.1 COSY-11 facility — general remarks 

The experiment described in this thesis has been performed at the COSY-11 facil- 
ity [13], an internal magnetic spectrometer installed at the cooler synchrotron COSY in 
Jiilich [14]. The COSY-1 1 detection system is schematically depicted in figure 3.1. De- 
tails of the functioning of all detector components and the method of measurement can 
be found in references [13, 15, 16]. Therefore, here the used experimental technique 
will be only briefly presented. 

The synchrotron accelerates protons and deuterons up to a momentum of 3.4 GeV/c. 
At the highest momentum a few 10 10 accelerated particles pass through the target ~ 10 6 
times per second. 

The hydrogen (H2) or deuteron (D2) cluster target (see figure 3.2) is installed in front of 
the dipole magnet. The positively charged products of the reaction are bent in the mag- 
netic field of the dipole and leave the vacuum chamber through thin exit foils, whereas 
the beam — due to the much larger momentum — remains on its orbit inside the ring. 
The charged ejectiles are detected in drift chambers (Dl, D2, D3) [13] and scintillator 
hodoscopes (SI, S2, S3) [13, 15]. Neutrons and gamma quanta are registered in the 
neutron detector (N). In order to separate neutrons and gamma quanta from charged 
particles veto detector (V) is used. An array of silicon pad detectors (Si spec ) is used 
for the registration of the spectator protons. Protons scattered under large angle are 
measured in another position sensitive silicon pad detector (Si mon ). 

The experiments performed at COSY-1 1 base on the measurement of four-momenta 
of the outgoing particles. Unregistered short lived mesons and hyperons are identified 
via the missing mass technique. 

For each charged particle, which gave signals in drift chambers, the momentum 
vector can be determined. First the trajectories of the particles are reconstructed [17], 
and then knowing the magnetic field of the dipole, the momentum vector is recon- 
structed. In case of two close tracks, the information about the energy loss from SI, S2, 
and S3 is used to inspect the efficiency of the track reconstruction. Particle's velocity 
determination is based on the time- of- flight measurement between SI (or S2) and S3 
detectors. Knowing the velocity and the momentum of the particle, its mass can be 
calculated, and hence the particle can be identified. After the particle identification the 
time of the reaction at the target is obtained from the known trajectory, velocity, and 
the time measured by the S 1 detector. The neutron detector delivers the information 
about the time at which the registered neutron or gamma quanta induced a hadronic or 
electromagnetic reaction. 
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Figure 3.1: Scheme of the COSY-11 detection system. Protons are registered in two drift 
chambers Dl, D2 and in the scintillator hodoscopes SI, S2, S3. An array of silicon pad detectors 
(Si spec ) is used for the registration of the spectator protons. Neutrons are registered in the 
neutron modular detector (N). In order to distinguish neutrons from charged particles a veto 
detector (V) is used. Deuterons with larger momentum are registered in deuteron chamber D3. 
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Figure 3.2: Schematic view of the cluster target. 



The time of the reaction combined with this information allows to calculate the 
time-of-fiight (TOF^) of the neutron (or gamma) between the target and the neutron 
detector, and — in case of neutrons — to determine the absolute value of the momentum 
(p) what can be expressed as: 



where m denotes the mass of the particle, / stands for the the distance between the target 
and the neutron detector and TOF N is the time-of-fiight of the particle. 

In order to calculate the four-momentum of the spectator proton — in case of quasi- 
free reactions with proton beam and deuteron target — its kinetic energy (T) is directly 
measured as the energy loss in the silicon detector (Si spec ). Knowing the proton kinetic 
energy one can calculate its momentum using the following relationship: 



where m denotes the proton mass, and T is the kinetic energy. In case of the measure- 
ments with deuteron beam and proton target the trajectory of the spectator proton may 
be reconstructed from signals from the drift chambers Dl and D2 and hence in this case 
its momentum can also be determined. 

To evaluate the luminosity, the elastically scattered protons are measured at the 
same time. With one proton detected in the drift chambers and scintillator hodoscopes 
and the other proton in the silicon detector Si mon resulting in the determination of the 
hit position, the elastically scattered protons can be well separated. 
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(3.1) 




(3.2) 
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3.2 Functioning of the neutron detector 

In this section I would like to emphasize the neutron detector, since the time calibration 
of this detector constitutes one of the main goals of this thesis. 

Previously, the neutron detector was built out of 12 detection units, with light guides 
and photomultipliers mounted on one side of the module. In order to improve the time 
resolution of the detector additional light guides and photomulipliers were installed, 
such that the light signals from scintillation layers are read out at both sides of the 
module. The present neutron detector consists of 24 modules, like the shown in figure 
3.3. Each module is built out of eleven plates of scintillator material with dimensions 
240 mm x 90 mm x 4 mm interlaced with eleven plates of lead with the same dimen- 
sions. The scintillators are read out at both edges of the module by light guides — 
made of plexiglass — whose shape changes from rectangular to cylindrycal, in order to 
accumulate the produced light on the circular photocathode of a photomultiplier. 




Figure 3.3: Schematic view of a single module of the neutron detector. LG and PM denote light 
guides and photomultipliers, respectively. Picture is taken from [18]. 




Figure 3.4: Photo of single modules of the neutron detector, at the stage when the scintilla- 
tor/lead structure was well visible. 

The neutron detector is positioned at a distance of 7 m behind the target in the con- 
figuration schematically depicted in figure 3.5. As can be deduced from figure 3.6 the 
maximum efficiency, for a given total thickness, for the registration of the neutron — 
in the kinetic energy range of interest for the COSY-1 1 experiments (~ 300 MeV - ~ 
700 MeV) — would be achieved for the homogeneous mixture of lead and scintillator. 
In order to optimize the efficiency and the cost of the detector the plate thickness has 



Experimental setup 



17 



been chosen to be 4 mm. This results in an efficiency which is only few per cent smaller 
than the maximum possible. 
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Figure 3.5: Schematic view of the neutron detector from above. The squares with numbers 
represent single detection modules. The z axis is defined by the beam direction. 
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Figure 3.6: Calculated efficiencies of the 20 cm thick neutron detector as a function of neutron 
kinetic energy for various layer thicknesses of the scintillator and lead plates. The solid line 
with triangles is for pure scintillator. Diamonds denote a layer thickness of 0.5 mm, closed 
circles of 2 mm, open circles of 5 mm, and crosses of 25 mm. The figure has been adapted 
from [19]. 
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The functioning of the detector was already confirmed in experiments carried out 
with a deuteron beam and hydrogen target. Figure 3.7 shows experimental distributions 
of the number of hits per individual detection unit, for neutral (left) and charged parti- 
cles (right). As expected the counting rate of modules in the middle part of the detector 
is much higher than these for the modules on the edges. In particular the smallest rate 
is observed in modules No. 5, 10, 15, and 20, because this row is partly out of the 
geometrical acceptance of the dipole yoke. The number of hits per segment from ex- 
periment is in perfect agreement with the corresponding spectra which were simulated 
using a GEANT-3 code (see figure 3.8). 
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Figure 3.7: Number of hits per individual module. Experimental spectra obtained for neutral 
particles (left) and charged particles (right) are shown. 
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Figure 3.8: Number of hits per individual module from Monte-Carlo studies obtained for 
neutral particles (left) and charged particles (right). The simulation was performed for the 
dp — ► ppn sp reaction taking into account reactions of the ejectiles with the dipole and all 
detector materials. 



4. Calibration of the neutron detector 



The installation of the neutral particle detector at the COSY-1 1 facility enables to study 
a plethora of new reaction channels. This detector is designed to deliver the time at 
which the registered neutron or gamma quantum induced a hadronic or electromagnetic 
reaction, respectively. This information combined with the time of the reaction at target 
place — deduced using other detectors — enables to calculate the energy of the detected 
neutron. In this section a method of time calibration will be demonstrated and results 
achieved by its application will be presented and discussed. Information about the 
deposited energy is not used in the data analysis because the smearing of the neutron 
energy determined in this manner is by more than order of magnitude larger than this 
established from the time-of-flight method. 



4.1 Time calibration of the neutron detector 



4.1.1 Time signals from a single detection unit 



As already discussed in section 3.2 the neutron detector at the COSY-1 1 facility is built 
out of 24 modules such as the one shown in the figure below. 

The time (T exp ) from a single module is calculated as an average time measured by 
the upper and lower photomultiplier. Namely: 



where Ttdc denote the time difference between the arrival of the photomultiplier and 
trigger signals to the Time-to-Digital-Converter (TDC). 




-.Up . rpd W 

TDC X TL) C 
2 



(4.1) 
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n 




Figure 4.1: Definition of variables x and L used in the text. The scheme of the detection module 
is the same as in figure 3.3. 



This can be expressed as: 



^TDC — treat ~l~ offset ^ -\- ^trigger-, 

CL 



(4.2) 



T^ c = t rea i + off set dw + 



L — x 
cl 



1 trigger 5 



(4.3) 



where L stands for the length of a single module, x denotes the distance between the 
upper edge of the active part of the detector and the point at which a neutron induced the 
hadronic reaction, t rea \ is the time at which the scintillator light was produced, T tr i gger 
represents the time at which the trigger signal arrives at the TDC converter, and cl 
denotes the velocity of the light signal propagation inside the scintillator plates. The 
parameters offset" p and offset^' denote the time of propagation of signals from the 
upper and lower edge of the scintillator to the TDC unit. 

Applying equations 4.1, 4.2, and 4.3 one can calculate a relation between T exp and 
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treal '■ 

offset"* + offset*** + ^ 

^ ^ = ^rea/ H 2 Ttrigger = treal "I" offset — Ttrigger (4.4) 

The value of "offset" comprises all delays due to the utilized electronic circuits, and 
it needs to be established for each segment. It is worth to note, that this time of the 
neutron detector signal is independent of the hit position, as it can be deduced from 
equation 4.4 and was proven experimentally, and depends on the difference between 
the time of light generation and the trigger only. 



4.1.2 Relative timing between modules 

Instead of determining the value of "offset" from equation 4.4 for each detection unit 
separately, the relative timing between modules will be first established and then the 
general time offset connecting the timing of all segments with the other detectors of the 
COSY-1 1 setup will be found. In order to establish relative time offsets for all single 
detection units, distributions of time differences between neighbouring modules were 
derived from experimental data. A time difference measured between two modules can 
be expressed as: 

Aij = T™ p - Tf p = tr l - t r / al + (offset; - offset;) , (4.5) 

where T t exp and Tj xp stand for the time registered by the i th and j th module, respec- 
tively. Example of A ;J spectra determined from a measurement carried out in January 
2003 with a hydrogen target and a deuteron beam accelerated to the momentum of pt> = 
3.204 GeV/c are presented in figure 4.2. The time differences were calculated assuming 
that all constants (offset) are equal to zero (see eq. 4.5). One can note that the peaks are 
shifted and additionally the distributions contain long tails. The tails reflect the veloc- 
ity distribution of the secondary particles. Corresponding spectra of time differences 
between the modules, shown in the figure 4.3, were generated using a GEANT-3 code. 
To produce these spectra the quasi-free dp — > ppn sp reaction has been simulated. The 
details are described in the appendix A. The values of the relative time offsets were 
determined using a dedicated program written in the Fortran 90 language [22, 23]. It 
adjusts values of offsets such that time differences obtained from experimental data and 
from simulation equals to each other for each pair of detection units. Furthermore, from 
the width of the spectra one can receive the information about the time resolution of a 
single module, which was found to be 0.4 ns [22]. 
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Figure 4.2: Distribution of the time difference between the 6 and the 11, the 1 th and the 12 , 
and the S th and the \3 th module of the neutron detector, as determined before the calibration. 




Figure 4.3: Simulated distribution of the time difference between the 6 and the 11 , the 7' 
and the \2 th , and the 8 th and the 13 th module of the neutron detector. 



nth 




Figure 4.4: Distribution of the time difference between the 6 and the 11, the 1 th and the 12 th , 
and the S th and the \3 th module of the neutron detector, as obtained after the calibration. 

Figure 4.4 presents experimental distributions of time differences between neigh- 
bouring modules as determined after the calibration. Examining figure 4.3 and 4.4 it is 
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evident that now the peaks are positioned precisely at the value expected from simula- 
tion. 

4.1.3 General time offset 

For the determination of the momentum of neutrons eg. in the analysis of reactions like 
dp — > ppny, the time between the reaction moment and the hit time in the neutron 
detector has to be determined for each event. The time of the reaction can be deduced 
from the time when the charged particle crosses the SI detector (see fig. 3.1) assuming 
that this particles trajectory and velocity can be reconstructed. To perform the calcu- 
lation of the time-of-flight between target and neutron counter a general time offset 
of the neutron detector with respect to the SI detector has to be established. For this 
purpose the quasi-free dp — > ppn sp reaction will be used. In this type of reactions the 
proton bound in a deuteron scatters elastically on a target proton, whereas the neutron 
considered as a spectator does not interact with the proton, but escapes untouched and 
hits the neutron detector. 

Data have been taken at a beam momentum of 3.204 GeV/c close to thereshold of 
the dp — > dpr\ process. Events corresponding to the dp — > ppn sp reaction have been 
identified by measuring the outgoing charged as well as neutral ejectiles. 



7 m 




Figure 4.5: Schematic view of the COSY-11 detection system. Only detectors used for the 
measurement of elastic scattered protons and spectator neutrons are shown. 

Fast protons are detected by means of the drift chambers (D) and scintillator ho- 
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doscopes (S1-S3). Protons scattered under large angles are measured in a position sen- 
sitive silicon detector (Si mon ). Neutrons are registered in the scintillator-lead sandwich 
detector (N). 

The time-of-flight between the target and the neutron detector (TOF neut ) is calcu- 
lated as a difference between the time of the module in the neutron detector which fired 
as the first one (t™ al ) and the time of the reaction {t r r eal ). 



T0F neut =f eal _ t real (4 6) 



Time of the reaction is obtained from backtracking of protons through the known mag- 
netic field and the time measured by the SI detector (T S1 ). Thus f r eal can be expressed 
as: 

t r r eal = tsi ~ TOF = T sl - offset 51 + T trigger - TOF, (4.7) 



where TOF denotes the time-of-flight between the SI counter and the target, and 
offset 51 denotes all delays of signals from the S 1 detector. Since time in the neutron 
detector reads: 

f n eal = T n - offset" + T trigger (4.8) 



we have: 



T0F neut = Tn + T0F _ T si _ ff set ^ + ff set si = T n + TOF - T S1 + offset (4.9) 



By offset G the general time offset of the neutron detector in respect to S 1 is denoted. 
In order to establish this global time offset, the time-of-flight spectrum derived from 
experimental data for the dp — > ppn sp reaction was compared with the corresponding 
distribution which was reconstructed from the signals simulated in the detectors (see 
figure 4.6 left). To arrive at the same statistic as was achieved in experiment, 8 • 10 
events for the dp — > ppn sp reaction were simulated using a GEANT-3 code. In the 
simulation the time resolution of a single module of o = 0.4 ns was taken into account. 
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Figure 4.6: (left) Time-of-flight distribution between the target and the neutron detector as 
obtained before the general time offset had been determined. The dashed histogram denotes 
the experimental data whereas the solid one depicts the result of the simulation, (right) Ex- 
perimental time-of-flight distribution (dashed histogram) calculated between the target and the 
neutron detector with general time offset equal to 13.6 ns compared to the Monte-Carlo sim- 
ulation (solid histogram). The latter was normalized in number of counts to the experimental 
data. 

The simulated spectrum was normalized so that the integrals of both distributions 
are equal. With a general time offset of 13.6 ns the experimental distribution corre- 
sponds to the simulated one as shown in figure 4.6 (right). The main cause of the 
smearing of the considered time distribution is the Fermi momentum of the nucleon 
inside the deuteron. The time resolution and dimensions of the detector are of minor 
importance. 



4.2 Momentum resolution 

The experimental resolution of the missing mass determination eg. in the analysis of 
reactions like pn — > pnx\' strongly rely on the accurate measurement of the momentum 
of the neutrons [9] . Therefore, the momentum resolution of the neutron detector has 
to be elaborated. After the neutron and gamma quanta are identified the momentum 
of neutrons is calculated from the time-of-flight between the target and the neutron 
detector according to the formula 3.1. 

Monte Carlo studies of the dp — > ppn sp reaction have been performed in order to 
establish the momentum resolution of the neutron detector. Figure 4.7 presents the dif- 
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ference between the generated neutron momentum (P gen ) and the reconstructed neutron 
momentum from signals simulated in the detectors (P re c)- 



AP = P, 



gen 



(4.10) 



The value of (P rec ) was calculated taking into account the time resolution of the neutron 
detector (a = 0.4 ns) as well as the time resolution of the SI counter (o = 0.25 ns). 




-1 -0.75 -0.5 -0.25 0.25 0.5 0.75 1 

A Q [GeV/c] 



Figure 4.7: Difference between generated (P gen ) and reconstructed (P rec ) neutron momenta for 
the dp — > ppn sp reaction simulated at deuteron beam momentum of 3.204 GeV/c. 

The distribution of AP was fitted by a Gaussian function resulting in a momentum 
resolution of o(P) = 0.14 GeV/c. Consequently the fractional momentum resolution 
g(P)/P for neutrons with momentum value of 1.6 GeV/c is equal to 8%. This fractional 
resolution changes with the momentum of the neutron (see Appendix B) and eg. for 
the neutrons produced at threshold for the pn — > pnr\ reaction it amounts to 3%. 



5. Analysis of the experimental data 



For the first time at the COSY-1 1 experiments signals from y-quanta were observed in 
the time-of-fiigth distribution (for the neutral particles) measured between the target 
and the neutral particle detector. A corresponding spectrum is shown in figure 5.1 (left). 
The data are from an experiment carried out using a deuteron target and a proton beam 
with a momentum of 2.075 GeV/c [9]. In addition to a broad distribution originating 
from neutrons, a sharp peak from y rays is seen at a value of about 25 ns. A Monte 
Carlo simulation performed for the pn — > pnr\ reaction, which is one of the possi- 
ble processes contributing to the neutron time-of-flight distribution, is shown in figure 
5.1 (right ) 




Figure 5.1: Time-of-flight distribution determined between the target and the neutron detector, 
(left) Experimental spectra, and (right) Monte Carlo simulation. Figure adapted from [20] 

In this chapter results of the analysis aiming for the identification of bremsstrahlung 
radiation in the data taken in January 2003 will be presented. The measurement was 
carried out with a deuteron beam and hydrogen cluster target at the cooler synchrotron 
COSY-Julich by means of the COSY-1 1 detection system, and its primordial aim was 
to investigate the dp — > dpr\ reaction close to the threshold [5]. The experiment was 
performed at four different deuteron beam momenta between pbeam = 3.165 and 3.204 
GeV/c. During the run with a beam momentum of Pbeam — 3.204 GeV/c, an additional 
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trigger with neutron detector — referred to T8 — was set up for the registration of 
charged ejectiles in coincidences with neutrons or gamma quanta. These conditions 
can be written symbolically as: 

T8 = Sl M=2 AN lu > 2 AS3^ =h 

which means that two signals in the SI detector, one signal in the S3 detector and, at 
least two signals in the neutron detector were demanded. 

Possible reactions with gamma quanta in the final state can be divided into two 
groups, viz: free dp — > dpy and dp — > 3 Hey reactions, and quasi-free dp — > dp sp y, 
dp — > ppjn sp , dp — > pnyp sp reactions. In case of quasi-free reactions one of the 
nucleons bound in a beam deuteron is treated as a spectator (p sp , n sp ) and does not take 
part in the reaction. 

As a first step of the data analysis events with simultaneous signals in any of the 
drift chambers and the neutron detector were selected. Figure 5.2 presents the time- 
of-flight distributions between the target and the neutron detector obtained assuming 
that in coincidence with a neutral particle also a proton (left) or deuteron (right) was 
identified based on signals from drift chambers and scintillator hodoscopes. The signals 
from y-quanta do not appear in the distributions which are predominantly due to quasi- 
free elastic dp — > ppn sp and dp — > dnn + reactions. 




Figure 5.2: Experimental time-of-flight distributions determined between the target and the 
neutron detector obtained under the condition that additionally one proton (left) or one deuteron 
(right) was registered in drift chambers. 

The charged ejectiles can be well identified as shown in the figure 5.3. Three clear 
peaks evidently visible in this figure correspond to the squared mass of pion, proton, and 
deuteron. The mass of the particle is calculated from its momentum — reconstructed 
from tracking back through magnetic field to the target point — and velocity determined 
from the time-of-flight measured between SI and S3 detector. 
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Figure 5.3: Distribution of the squared mass of charged particles originating from the 
deuteron-proton reaction performed at a beam momentum of 3.204 GeV/c. 



5.1 Identification of the dp —> dpy reaction 



In order to identify the dp — > dpy reaction events with two tracks in drift chambers and 
a simultaneous signal in a neutron detector have been selected. In figure 5.4 squared 
mass of one particle is plotted versus squared mass of the other registered particle. Base 
on this figure measured reactions can be grouped according to the type of ejectiles. 
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Figure 5.4: Scatter plot of invariant masses determined for events with two charged particles 
measured in coincidence. 



Thus reactions with two protons, proton and pion, proton and deuteron, and pion 
and deuteron can be cleary separated. Figure 5.5 shows the experimental distribution 
of the time-of-fiight between the target and neutral-particle detector with the require- 
ment that two charged particles were registered and that one of them was identified as a 
proton and the other as a deuteron. In this case due to the baryon number conservation, 
there is only one possible source of a signal in a neutron detector, namely a gamma 
quantum. In fact a clear peak around the time of 24.5 ns is visible, and this is just the 
value corresponding to the time-of-flight of light on a distance of 7m. The gamma 
quanta may originate from the bremsstrahlung reaction (dp — > dpy) or from the decay 
of mesons produced eg. via dp — > dpn° — > dpyy or dp — > dpn°n° — > dp4y reac- 
tions. It is possible to distiguish between these hypothesis calculating the missing mass 
produced in the dp — > dpX reaction. 
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Figure 5.5: Time-of-flight distribution determined between the target and the neutron detector 
obtained under the assumption that additionally to a signal in the neutron detector, one proton 
and one deuteron were identified from the signals in the drift chambers. 

Knowing the four momenta of a proton and a deuteron in the initial and final state, 
and employing the principle of momentum and energy conservation one can calculate 
the squared mass of the unmeasured particle or group of particles: 

m\ = E 2 X - p 2 = (E b + E t -E d -E p ) 2 - (p b + p t -p d -p p ) 2 (5.1) 

where, 

Eb,pb is the energy and momentum of deuteron beam, 
E t ,p t is the energy and momentum of proton target, 
Ed,Pd is the energy and momentum of outgoing deuteron, and 
Ep,p p is the energy and momentum of outgoing proton. 

Figure 5.6 shows the squared missing mass distribution as obtained for the 
dp — > dpX reaction. A significant peak around MeV 2 /c 4 — the squared mass of 
the gamma quantum — constitutes an evidence for events associated to the deuteron- 
proton bremsstrahlung (dpy). In addition a broad structure at higher masses originating 
from two pions emitted from the reaction dp — > dpn n° is visible. 
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Figure 5.6: Squared missing mass distribution for the dp — > dpX reaction. 



One of the most important information is the energy spectrum of the gamma quanta. 
The energy of gamma quanta can be measured with a very poor accuracy, however it 
can be calculated using a simple relation, after selecting events corresponding to the 
dp — > dpy reaction, namely: 

E y = E b + E t -E d -E p (5.2) 

where, 

E[, is the energy of deuteron beam, 

E t is the energy of proton target, 

E d is the energy of outgoing deuteron, and 

E p is the energy of outgoing proton. 

Figure 5.7 shows the determined spectrum. One sees that the registered gamma 
quanta populate predominantly the energy range between 0.8 and 0.9 GeV, which is 
partially due to the acceptance of the COSY-1 1 system which decreases rapidly with 
increasing kinetic energy shared by the ejectiles. The detailed conclusions concerning 
the real energy distribution of the produced gammas will require careful acceptance 
corrections which will be performed in the near future. At present we can consider the 
observed distribution as an evidence that the produced quanta are high energetic. 
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Figure 5.7: Experimental distribution of gamma energy for the dp — > dpy reaction. The 
spectrum is not corrected for the detection acceptance. 
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6. Summary and perspectives 



The experiment which was described in this thesis has been performed at the cooler 
synchrotron COSY in the Research Center Jiilich by means of the COSY-1 1 detection 
system. 

The results of the identification of bremsstrahlung radiation in data taken with a pro- 
ton target and a deuteron beam have been presented and discussed. For the first time 
— using COSY-1 1 facility — events associated with the dp — > dpy reaction have 
been observed. For the quantitative determination of the total cross section of the 
dp — > dpy reaction the luminosity and detection acceptance remains to be established. 
There are also plans to analyse the data in view of bremsstrahlung radiation in a quasi- 
free np — > npy reaction. 

The second point of this work was the time calibration of the neutron detector. As 
shown in the chapter 4, the general time offset of neutron detector with respect to the 
SI detector was found to be 13.6 ns. The resolution of the neutron momentum de- 
termination by means of the neutral particle counter — the crucial factor in neutron 
momentum determination — was found to be 0.14 GeV/c at a neutron momentum of 
1.6 GeV/c, the dependence of the fractional momentum resolution as a function of a 
neutron momentum is presented in Appendix B. 



The installation of the neutron detector enables not only to study the isospin depen- 
dence of the meson production [9] and bremsstrahlung radiation, but also to investigate 
the production of the resonance + in the elementary proton-proton interaction. A 
signature of the + production may be the presence of a 1.54 GeV/c 2 peak in the n K + 
invariant mass distribution for the pp — > nK + L + reaction. The data analysis aiming 
for the determination of the n K + invariant mass distribution has just started. In order 
to determine the acceptance of the COSY-1 1 detection system for the pp — > nK + L + 
reaction we have simulated the response of the detectors for 6 ■ 10 6 events generated 
in the target. The solid histogram in figure 6. l(left) illustrates the distribution of the 
invariant mass nK + for the generated events while the dashed line depicts the spectrum 
which was reconstructed from the signals simulated in the detectors. 
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Figure 6.1: (left) Phase space distribution of the invariant mass of the nK + system of the 
pp — > nK + L + reaction (solid line) and its convolution with the COSY-11 acceptance (dashed 



line), (right) The acceptance of the COSY-11 detection setup for the pp 
tion [21]. 



nK + V 



reac- 



The ratio of the obtained invariant mass distributions results in the differential ac- 
ceptance of COSY-11 facility for detecting the pp — > nK + Z + reaction as shown in 
figure 6.1 (right). 

The total acceptance of COSY-1 1 detection system for the pp — > nK + L + reaction 
measured at the beam momentum of 3.257 GeV/c is equal to 10~ 4 [20]. 

Additionally, the simulation of the invariant mass distribution of the process 
pp — > E + + — > nK + Z + has been performed, taking into account the width of + 
equal to 5 MeV [24] . The results of the Monte-Carlo calculation are shown in figure 
6.2(left). 



The expected signal from the pp — > E + ® + — > n K + L+ reaction together with the 
background originating from the direct pp — > nK + Z + reaction is presented in figure 
6.2(right). 

Here it is assumed arbitrarily that the total cross section for the pp — > E + + reaction 
is ten times smaller than the one for pp — > nK + Z + . If the performed appraisals are 
realistic, one should observed a clear signal in the experiment originating in the © + 
production as it is noticeable in the right panel of figure 6.2. 
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Figure 6.2: (left) Invariant mass distribution of the nK + system for the 
pp — > Z + + — > reaction, (right) The expected full invariant mass distribution of 

the nK + system [21]. 
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A. Kinematics of the dp — > ppn sp 
reaction 



The quasi-free reaction dp — > ppn sp has been used to determine (i) the relative timing 
between modules, (ii) the general time offset of the neutron detector, and (iii) to estab- 
lish the momentum resolution of this detector. Due to the Fermi motion of the nucleons 
bound in the deuteron the simulation of the quasi-free reaction proceeds in following 
steps : 

The components of the Fermi momentum of a nucleon in Cartesian coordinate system 
are equal to: 

Pfx = Pf ■ sin0cos(|) 
Pf y = pf • sin0sin(|) 
Pf z = Pf ■ cos 9 

As a first step the value of pf (absolute Fermi momentum) is generated according 
to the momentum distribution of nucleon inside the deuteron derived from the PARIS 
potential model [18, 25] . Next the azimuthal angle (]) and cosine of polar angle © (cos ©) 
which define the momentum direction, are generated assuming an uniform distribution. 
Further the nucleon Fermi momentum inside the deuteron is related to the nucleon 
momentum in laboratory frame by Lorentz transformation [26]: 

p/ ab =p/+ pvky/(y+ m-p f + E f ) 

Ef b = y(E f + %.p f ) 
where is the velocity of the deuteron in the laboratory frame and y is equal to: 

After the momenta of the nucleons inside the deuteron are converted into the LAB 
system, the proton from the deuteron scatters elastically on a proton target, whereas 
the neutron does not take part in the reaction, but with momentum possessed at the 
time of the reaction remains untouched. In order to simulate the proton-proton elastic 
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scattering now we calculate the proton momentum in the proton-proton system using 
the transformation equation: 



\p cm \ = ^(M 2 -4m 2 p )/4 
where M is the total mass of the colliding protons: 



M=J(Jp/ ab2 +ml + m p r-p f lah2 



Once more cos & cm and § cm , which define the momentum direction of protons after 
the scattering in the proton-proton system are generated assuming an uniform distribu- 
tion. Thus, the components of the proton momentun after scattering are equal to: 

p cm' = \ p Sm\ s i n Qcm cos ^cm 

pf = |p cw |sin0 cw sin(|) cm 
pf = \p cm \cosQ cm 

As a last step the proton momentum after scattering in the proton-proton system is 
related with momentum in laboratory frame by Lorentz transformation: 

pfab' — pcm' _|_ (3?ray:™ (y™ j (*fm _|_ 1 ) pcm . p 7m' _|_ gem' j ^ 

where 



(3™ 



P f 



p'f +m 2 p + m p 



and 



-> 2 
Jab 



pf b +m 2 p + m p 



M 



B. Resolution of the measurement of 
the neutron momentum 



The momentum resolution of the neutron detector is a crucial factor in the neutron mo- 
mentum determination and as shown in figure 6.3 it strongly depends on the momentum 
of the neutron. 
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Figure 6.3: Fractional momentum resolution of the neutron detector as a function of neutron 
momentum. 



The momentum of the neutron is calculated from the time-of-flight between the 
target and the neutron detector and can be expressed as: 
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where m denotes mass of the particle, I stands for the distance between the target and 
the neutron detector and t N is the time-of-flight of the particle. A fractional momentum 
resolution is given by the equation: 

dp _ 

p p 

where o r accounts for both the time resolution of the neutron and S 1 detectors and 
can be written as: 



<*t = \/ 



SI 



The result presented in figure 6.3 was obtained assuming that o„ = 0.4 ns, 
Gsi = 0.25 ns, and 1 = 1 m. 



C. Dalitz plot distribution for the dpy 
system at = 3.37 GeV 



As a first step for establishing of the double differential acceptance of the COSY-11 
facility for the measurement of the dp — > dpy reaction a Dalitz plot distribution for the 
dp — > dpy reaction was calculated. Figure 6.4 presents the result of calculations as- 
suming that the phase space volume is homogeneously populated. Modifications of this 
distribution due to the COSY-1 1 detection setup acceptance remain to be determined. 




Figure 6.4: Dalitz plot distribution as generated for the dp — ► dpy reaction. 

Squared of the invariant masses of the py and dy system can be calculated according 
to the below equations: 

Spy = (E p +Ey) 2 ~ (p p + Py) 2 
43 
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Sdy = (Ed + E y ) 2 - (p d + p y ) 2 
It is interesting to note that minimum values for the S p y and S^y are equal to: 

py P 

c<min is! 

^dy ~ M d 
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